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Abstract. We study the power spectra of the variability of seven intermediate polars containing magnetized 
asynchronous accreting white dwarfs, XSS J00564+4548, ICR J00234+6141, DO Dra, V1223 Sgr, IGR J15094- 
6649, IGR J16500-3307 and IGR J17195-4100, in the optical band and demonstrate that their variability can be 
well described by a model based on fluctuations propagating in a truncated accretion disk. The power spectra 
have breaks at Fourier frequencies, which we associate with the Keplerian frequency of the disk at the boundary 
of the white dwarfs' magnetospheres. We propose that the properties of the optical power spectra can be used to 
deduce the geometry of the inner parts of the accretion disk, in particular: 1) truncation radii of the magnetically 
disrupted accretion disks in intermediate polars, 2) the truncation radii of the accretion disk in quiescent states 
of dwarf novae. 
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1. Introduction 

Aperiodic variations in the intensity of many astrophysi- 
cal sources carry a lot of information about physical pro- 
cesses in and around them, which is complementary to 
and completely independent of all other types of astro- 
physical information. It was noticed already long ago that 
accreting sources typically demonstrate flux variability on 
a wide range of time scales, from milliseconds and seconds 



was recognized that photons, which contain this flux vari- 
ability, originate in the innermost regions of the accretion 
flow. While for X-ray accreting sources it was quite clear 
from the beginning (simply from energy contained in the 
variable part of the X-ray radiation), for accreting white 
dwarfs it was shown with the help of eclipse mapping (e.g. 
Bruch||1992[|1996||Baptista k Bortoletto||2004[ ) . 



However, until recently the origin of the these vari- 
ations remained unclear. This variability possesses a set 
for galactic compact objects (e.g. Linnell||1950 for accret- of properties that was not easy to explain in early mod- 
ing white dwarfs, |Rappaport, Doxsey, k Zaumen|[T971[ els. For example, the shot-noise model (see e.g. Terrell 
Oda et al. 1974 for X-ray binaries) to months and years 1972 ) proposed that the variability originates as a result 



for accreting supcrmassivc black holes in active galactic 
nuclei (e.g. Zaitseva k Lyutyi]|1969[ ). Aperiodic flux vari- 
ations of X-ray emitting compact sources were discov- 
ered since the 1970s and have been extensively studied 
since then. Typically these aperiodic flux variations have 
a power spectrum (square of Fourier transform of the light 
curve) that can be described as a power law with a slope 
P oc y- 1 ----i-5 ) sometimes with flattening at low Fourier 
frequencies and steepening at high Fourier frequencies. It 



of additive summation of randomly occurring exponential 
"shots". This model explained the shape of power spec- 
tra of some X-ray binaries (e.g. Terrell|1972 Kraicheva et 



al. 1999), but failed to reproduce very essential observa- 



tional properties: 1) the extremely wide range of variabil- 
ity time scales, on which the variability has a self-similar 



power spectrum in some cases (Churazov, Gilfanov, k 



Revnivtsev 2001), 2) the linear dependence of the am- 



plitude of flux variability on the time average flux of the 
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source (Uttley & McHardy 2001), 3) the log- normal dis- 
tribution of instantaneous values of the flux and 4) the 
behavior of the frequency dependent phase lags between 



lightcurves in different energy bands ( 


Miyamoto et al. 


1988 


Kotov, Churazov, & Gilfanov||2001 


)■ 



Now it is widely accepted that the most appropri- 
ate model of the flux variability of accreting sources is 



the model of propagating fluctuations ( 


Lyubarskii 


1997 


Churazov, Gilfanov, & Revnivtsev 


2001 


). In this model 



the variability of observed flux originates as a transfor- 
mation of variability of the instantaneous mass accretion 
rate in the innermost regions of the accretion flow (where 
the majority of the emission is created) into a variabil- 
ity of the emergent luminosity. In turn, the variability of 
the mass accretion rate in the innermost region is a result 
of propagation (in the flow going inward) of variability 
created at all outer radii of the extended accretion disk. 
All inner regions of the disk create additional modula- 
tions to the instantaneous mass accretion rate on top of 
already existing longer time scale variability, transported 
to this radius from outer parts of the disk. The combined 
action of all radii of the disk creates power spectra with 



relatively shallow slopes (—1 . . . —1.5) (see e.g. 


Lyubarskii 


1997 


Churazov, Gilfanov, & Revnivtsev 2001 


). This rel- 



atively simple construction allows us to explain all the 
major observational characteristics of flux variability (see 
more details in |Churazov, Gilfanov, fc Revnivtsev 2001 



Kotov, Churazov fc Gilfanov| [2001] |Uttley fc McHardy| 
2001[|Revnivtsev||2008[ ). ~ 

In the framework of this model one can make very spe- 
cific predictions for the shape of the power spectra of flux 
variability of accretion flows. For example, the accretion 
disk with an inner boundary should have a break in the 
power density spectrum of its time variability, approxi- 
mately at the Fourier frequency, which corresponds to the 
typical time scales of variability introduced into the flow 
at the inner radius of the disk. 

Accretion-powered X-ray pulsars and asynchronous 
magnetic white dwarfs (intermediate polars) have mag- 
netic fields strong enough to disrupt the inner parts of the 
accretion disks. Thus the the fastest variability timescales 
associated with the innermost regions of the disk should 
be absent or reduced in their power spectra. 

This assumption was recently explored by Revnivtsev 



et al. ( 2009 ) in the X-ray energy domain, who analyzed a 
large sample of magnetized accretion- powered objects (X- 
ray pulsars and an intermediate polar) and showed that 
the variabilities of all these sources indeed have a break in 
their power spectra at the anticipated Fourier frequencies. 
The frequency of these breaks corresponds to the the inner 
disk Keplerian frequency at the boundary of the central 
object's magnetosphere. The conclusion was strongly sup- 
ported because the frequency of the break varied during 
the giant outbursts of accretion-powered X-ray pulsars, in 
accordance with the simple prediction of how the radius 
of the dipole magnetosphere of the central object (and 
consequently the disk Keplerian frequency at this bound- 
ary) should vary with the mass accretion rate (X-ray lu- 



minosity). Among other things, this study opens a new 
possibility to measure the properties of magnetosphcres 
of compact objects via the study of their flux variabil- 
ity. In a lot of cases this approach is much easier than 
trying to measure the cyclotron absorption lines in spec- 
tra of pulsars or polarization of photons originated in the 
magnetized accretion curtains in intermediate polars. 

We present results of our pilot study of aperiodic vari- 
ability of some asynchronous magnetized accreting white 
dwarfs (intermediate polars) with accretion disks (this 
is an essential ingredient of binary systems in our sam- 
ple, because sometimes asynchronous magnetized accret- 
ing white dwarfs might be stream-fed rather then disk-fed 
, see e.g. |Hellier fc Beardmore||2002[ ) in the optical bands 
(B,V,g',R) and show the applicability of the method. 

Observational data for our study was collected with 
the help of the Russian- Turkish telescope RTT150, the 
Southern African Large Telescope (SALT), 0.76m and lm 
telescopes of South African Astronomical Observatory. 

2. Observational data 

Power spectra of variability of optical light curves of 
accreting intermediate polars were obtained via Lomb- 



Scargle periodograms (Lomb 1976 Scargle 1982), aver 



aged over all available uninterrupted datasets. As the typ- 
ical length of datasets is approximately several hours, the 
lowest Fourier frequency probed by our power spectra were 
~ 10~ 4 Hz, while the highest frequencies depend on the 
time resolution used. In Table 1 we present the list of 
sources, used in our analysis. Optical magnitudes of bi- 
naries are approximate, because of significant long term 
flickering, present in all systems. Photometry of SALT is 
subject to additional uncertainty because of the method 
of observations (see below). 

2.1. RTT150 

Optical observations of the intermediate polars 
XSS J00564+4548, IGR J00234+6141 and DO Dra 
(=YY Dra) were carried out with 1.5-m Russian- Turkish 
Telescope (RTT150) at TUBTTAK National Observatory 
(TUG), Bakirlitcpc mountain, Turkey. The objects were 
observed as a part of a long-term program of study 
of aperiodic optical variability of cataclysmic variables 
(CVs). 

We used a low-readout noise back- illuminated 2x2 K 
Andor DW436 CCD in combination with B, g' , V, or R fil- 
ters, mounted in F/7.7 Cassegrain focus of the telescope. 
The measurements with > 1 s time resolution were ob- 
tained from the series of direct images taken with only a 
small CCD subframe near the object. The bias was sub- 
tracted and flat-field correction was applied for all im- 
ages. The subframe included also at least one bright non- 
variable reference star, used for differential photometry. 
The use of differential photometry allowed us to diminish 
the influence of atmospheric turbulence on the resulting 
measurements, which is extremely important for aperiodic 
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Table 1. Optical observations of accreting white dwarfs used for construction of power spectra. 



Target 


Instr. 


Dates 


Mag 


Time res., sec 


Exp., ksec 


Spin Period, s 


XSS J00564+4548 


RTT150 


2005.10.19-23,26,28 


14.9(V) 


17-128 


71.6 


465 






2005.12.05,06 


14.6(R) 








IGR J00234+6141 


RTT150 


2005.09.10-11,12.01 


16.7(R) 


60-70 


37.6 


561 






2005.12.10-11,2007.09.03 


16.7(R) 








DO Dra 


RTT150 


2007.04.24,27,28 


14.7(R) 


1 


48 


529 






2008.05.24-26,06.03 


14.8( S ') 


1-8 


39 




V1223 Sgr 


SALT 


2008.05.13,2008.06.12 


~13.3(B) 


0.1 


2.9 


745 


IGR J15094-6649 


SAAO, lm 


2008.04.26-29 


14.6(V) 


7-8 


14.9 


809 


IGR J16500-3307 


SAAO, lm 


2008.04.27-28,05.04 


15.9(V) 


8-10 


18.4 


598* 




SAAO, 0.76m 


2008.08.03-05 


15.9(V) 








IGR J17195-4100 


SAAO, 0.76m 


2008.08.5-10 


15.2(V) 


8 


32.3 


1139* 



* — here we use the optical pulsation period, which also may be the beat frequency instead of the true spin frequency of 
the white dwarf 



variability studies. Data reduction and photometric mea- 
surements were made with IRAjQand our own software. 

The measurements with « 1 s time resolution were 
made with the setup where the narrow strip of the CCD 
was binned in its width into only one pixel, so that only 
one line of the data was read out of the CCD in each 
exposure. That was done to diminish the CCD read-out 
time. This setup also allowed us to make the differential 
photometric measurements, since the bright non-variable 
bright reference star was also observed in another part of 
this strip. We checked that these differential photomet- 
ric measurements effectively eliminated the influence of 
atmospheric turbulence even at these higher frequencies. 
The detailed description of this setup and appropriate cal- 
ibration measurements can be found in Burenin et al. (in 
prep.). 

2.2. SAAO 

Light curves of three intermediate polars IGR J15094- 
6649, IGR J16500-3307 and IGR J17195-4100 were ob- 
tained with the SAAO 1-m and 0.76-m telescopes and the 
UCT CCD (this is a frame-transfer CCD, implying that 
there is no dead time between exposures). Detailed de- 



scription of these data is presented in Pretorius (2009) 



2.3. SALT 

High time-resolution observations of V1223 Sgr were taken 
at the prime focus of the newly available 10-m class 



Southern African Large Telescope (SALT; jBuckley, Swart 
fc Meiring||2006| |Q'Donoghue et al.|[2006| during its per- 
formance verification phase on 2008 May 13 and June 
12. Observations were done with the imaging camera 



(SALTICAM; |0'Donoghue et al.|2003[[2004l ) in 'slot mode' 
( |0'Donoghue et al.|2006[ ). The SALTICAM is a CCD mo- 
saic of two 2048x4102 pixels CCDs with full image size of 
about 9.6x9.6 arcmin, but the nominal science field is 8 
arcmin. In 'slot mode', a mask with a slot of the width of 
20 arcsec (144 unbinned rows) is projecting on the CCDs. 
Only the central 72 rows receive all the light. The slot ex- 
tends across the full 8 arcmin in the horizontal direction. 
After each exposure, only 144 rows are transferred across 
the frame-transfer boundary, which shortens the vertical 
clocking overhead to 14 ms. The instrument thus provided 
a sequence of images 8 arcmin long, but only 20 arcsec 
wide. These images were split into four distinct sections, 
corresponding to the four readout amplifiers of the two 
CCDs. 

Total tracking time was ^2700 sec on May 13 and 
~1300 sec on June 12. The observations were made with 
the B filter. In order to facilitate the shortest possible 
exposures, 6x6 binning, or 1.2 arcsec/pixel was used. The 
exposure time for each image was 112 ms (the time to 
move each image behind the mask, during which image 
smearing takes place, is 14 ms). The observed data were 
all referenced to the nearby star of similar magnitude, 63 
arcsec to the east and 7 arcsec south of the target star. The 
comparison star and the target star were imaged within 
the same amplifier. No flat-field calibration frames were 
obtained (the calibration system was not available at the 
time of the observations) and the frames were simply bias- 
subtracted with overscanned pixels from each row. 

The data reduction of "slot mode" was done with the 
the PySALT user package. This tool is primarily written in 
python/PyRAF with some additional IRAF cod^ 

Analysis of the lightcurves obtained with SALTICAM 
in the described mode revealed problems in power spec- 



1 http://tucana.tuc.noao.edu/ 



See 



http:/ / www.salt.ac.za/science-support/salt-data- 



reduction/pysalt-users-package/ for more information. 
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tra of celestial sources at Fourier frequencies higher than 
^lHz, while at lower frequencies the results remained 
reliable. Therefore in our subsequent analysis we have 
not considered power spectra of the considered source at 
Fourier frequency higher than 0.7 Hz. 

Power spectra of V1223 Sgr obtained in different obser- 
vations were virtually indistinguishable from each other, 
therefore we averaged them. 

3. Discussion 

3.1. Power spectra of optical emission of CVs 

Optical emission of white dwarfs^] consists of two parts: 
one is provided by the internal energy dissipation inside 
the disk and another is a reprocessed EUV or X-ray emis- 
sion of white dwarf polar caps. 

In the framework of the model of propagating fluctua- 
tions we predict that the variability of the optical light in 
intermediate polars will be similar to that of the X-rays. 
There are two main scenarios: 

1 . If the optical emission of the system originates from the 
reprocessed X-rays, the X-ray and optical variability 
should be similar by construction. 

2. If the optical light is due to internal heating of the ex- 
tended accretion disk, we should still expect to observe 
similar power spectra of optical and X-ray light curves. 
That is because in both cases the flux variability is 
caused by the same fluctuations of the mass accretion 
rate. These fluctuations first affect the optical emis- 
sion of the disk and after that, when this same matter 
falls down onto the polar caps of the white dwarf, it 
changes the X-ray emission in the accretion column. 

In both cases power spectra of X-ray and optical variabil- 
ity of the CVs should be similar, and we can study the 
shape of the power spectrum of the X-ray light curve via 
observations of the optical variability of these CVs. 

The majority of white dwarfs are faint in X-rays (typ- 
ical fluxes are less than on the order of few photons/1000 
sq.cm./sec), therefore it is very difficult to observe them in 
X-rays. However, they can be more easily observed in the 
optical band. Therefore the study of optical variability in 
accreting white dwarfs can provide important diagnostics 
of the geometry of the inner parts of the accretion flow. 
In particular, we can immediately apply the diagnostics 

— to make estimates of the magnetic moment of white 
dwarfs. In the simplest model of dipolar magneto- 
sphere of the white dwarf, the dependence of the in- 
ner disk radius on the mass accretion rate in the disk 
is determined mainly by the white dwarf magnetic 
moment. Therefore measurements of power spectra of 

3 We emphasize that we consider here only magnetic cata- 
clysmic variables with accretion disks which might not always 
be the case even for asynchronous magnetic CVs, taking into 
account the existence of stream fed systems like V2400 Oph, 
see, e.g. Hellier & Beardmorei 2002 




Fig. 1. Scheme of accretion flow in magnetized accreting 
white dwarfs 



white dwarf variability in different accretion states will 
allow us to make estimates of the magnetic field. The 
proposed approach might be more sensitive than the 
search for cyclotron emission lines in optical-NIR spec- 
tra of objects. 

to make estimates of holes in the innermost regions 
of the accretion disk due to evaporation, predicted in 
some models of the quiescent accretion on white dwarfs 
(e.g |Meyer fc Meyer-Hofmeister||1994[ ) . 



3.2. Observed power spectra 

The observations of one of the brightest intermediate po- 
lars, V1223 Sgr, agree well with the considerations pre- 
sented above. In Fig. [2] we present the power spectra of 
V1223 Sgr in optical and in X-rays (from Revnivtsev et 
al. 2009). The similarity of these power spectra is strik- 



ing. We fitted these power density spectra with a simple 
analytical function P(f) cx /^(l + [///o] 4 ) _1/4 (essen- 
tial property whose is that it has a slope P oc / _1 at 
frequencies f < fo and a slope P oc f~ 2 at higher Fourier 
frequencies). The values of the break frequencies in opti- 
cal and in X-rays are: /break, opt = (2.1 ± 0.5) x 10 -2 Hz, 
/break, x-ray = (3.36 ± 0.3) x 10 -2 Hz (in order to deter- 
mine the /break values and their confidence intervals we 
used the x 2 minimization technique). The break frequen- 
cies are compatible at the 2<r level. There is marginal ev- 
idence for a higher break-frequency in X-rays, but it may 
be explained, for example, by a long-term difference in the 
accretion-rate of the source, since the optical and X-ray 
data were obtained in different epochs. 
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Fig. 2. Power spectra of V1223 Sgr in X-ray and in optical 
spectral bands. Dashed curves denote the simple analytic 
model described in the text ( P oc + [///o] 4 ) _1/4 ) 

fitted to the X-ray (black curve) and optical (red curve) 
power spectra. 

Assuming that the break frequency corresponds to the 
Keplerian frequency at the boundary of the white dwarf 
magnetosphere ( Revnivtsev et al.||2009 ), we can estimate 
its radius as 

r m ~ 1.6 x 10 9 cm 



or ~ 2.8i?wD (here we adopted the mass of the white 
dwarf in this system is 



WD 



= O.95M , Suleimanov,| 



Revnivtsev, & Ritter 20051. However, we should keep in 
mind that due to a lack of exact understanding of the for- 
mation of the break that this estimate is only accurate to 
within a factor of a few. 

In Fig. [3] we present power spectra of optical variability 
of seven objects — XSS J00564+4548, IGR J00234+6141, 
V1223 Sgr, DO Dra, IGR J15094-6649, IGR J16500-3307 
and IGR J17195-4100, measured with data described in 
the Sect. [2] The frequency axis of the power spectra was 
multiplied by the spin period of the white dwarf magneto- 
spheres. In the power spectra of a majority of these objects 
we removed the peak (simply by removing measurements 
of a power at Fourier frequencies near the inverse of WD 
rotational period), which was caused by spin modulated 



Table 2. Parameters of the break frequencies of the power 
spectra of the studied sources. 



Source 


/break x Period 


XSS 00564+4548 


0.7 ±0.2 


IGR J00234+6141 


1.0 ±0.3 


DO Dra 


0.5 ±0.1 


V1223 Sgr 


16 ±4 


IGR J15094-6649 


1.3 ±0.3 


IGR J16500-3307 


1.5 ±0.2 


IGR J17195-4100 


2.0 ±0.3 



variations of the optical brightness of the sources. Fainter 
(and possibly wider) QPO peaks may be present also in 
power spectra of objects without us detecting them. If 
such QPO peaks reside near the break frequency, they 
could somehow contribute to the observed sharpness of the 
break in some power spectra in Figj3| However, the dis- 
cussion of these possible QPO peaks is beyond the scope 
of this paper. 

In the framework of our model we predict that on this 
plot all accreting magnetic white dwarfs, in which the in- 
ner parts of the accretion disks are in corotation with 
white dwarf magnetospheres, will have breaks in the power 
spectra at a value around unity, i.e. around the WD spin 
frequency (see e.g. Figj3] and Table 2). If the accretion 
disk continues significantly below the corotation radius, 
the break in the power spectrum should be observed at a 
higher Fourier frequency. 

For comparison, we also show in Fig |3] the power spec- 
tra of an accreting X-ray pulsar, A0535±26, during its low 
and high luminosity state, i.e. during the periods of its 



high and low mass accretion rates (from Revnivtsev et al. 
2009). One can see that the power spectrum of the X-ray 



emission of A0535±26 with a high mass accretion rate con- 
tinues to higher Fourier frequencies (/ x Spin Period > 1), 
while variability at lower frequencies (fx Spin Period & 1) 
remained at almost the same level. 

The difference at high frequencies can be explained as 
an addition of an extra noise component, which is presum- 
ably generated in the annulus of accretion disk between 
the radii corresponding to the magnetosphere size in the 
high accretion rate (small radius) and low accretion rate 
(large radius) . This ring and therefore the associated vari- 
ability were absent in the state with low accretion rate. 

We can see that the behavior of accreting WDs is qual- 
itatively very similar to that of accreting neutron stars. 
Power spectra of all IPs except for VI 223 Sgr indicate 
that their white dwarfs are close to corotation with the 
inner edges of their accretion disks, while in V1223 Sgr 
the accretion disk should be truncated at a radius smaller 
than that of corotation. Simple models predict a spin-up 
of the central object in this case. 

This prediction does not agree with the detection of 
spin down, reported for this system approximately 20 
years ago by |van Amerongen, A ugustcij n, fc van Para dijs 
(1987), but this discrepancy may be explainable. For ex- 
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Fig. 3. Power spectra of optical variability of a sample 
of intermediate polars. Power spectra of DO Dra, XSS 
J00564+4548, IGR J00234+6141, ICR J17195-4100, and 
IGR J15094-6649 were scaled down for clarity by fac- 
tors, written on the figure. For comparison we present 
two power spectra of accreting X-ray pulsar A0535+26, 
collected during episodes of high and low mass accretion 
rates (thick solid histograms). See text for the details. 

ample, the source may have started to spin up since that 
time. Therefore additional observations are needed to in- 
vestigate this. Also, there may be a more complicated in- 
teraction of the accretion disk with the white dwarf mag- 
netosphere, in which the angular momentum of the white 
dwarf is diminishing while the accretion disk extends sig- 
nificantly below the corotation radius (corotation radius 
in this case is ~ 16i?wD)- As an argument in favor of 
this scenario we can mention that there are many of cases 
when actively accreting magnetic neutron stars do have a 
spin-down, while in a simple picture of accretion the mat- 
ter settling to the surface of the star should always speed 
up the rotation (if the disk is prograde). 

4. Summary and perspectives 

We demonstrated that high quality power spectra of opti- 
cal and X-ray variability of the intermediate polar, V1223 
Sgr (accreting white dwarf with a magnetosphere that dis- 
rupts the accretion disk at relatively large distances from 
white dwarf surface), are very similar to each other. Both 
power spectra have a very distinct break at a Fourier fre- 
quency which we associate with the Keplerian frequency 
at the inner radii of the truncated accretion disk. At higher 
Fourier frequencies the power spectra steepen. 

The similarity of the V1223 Sgr power spectra in the 
X-ray and the optical is expected in both scenarios of the 
origin of optical emission — if optical emission is domi- 
nated by the reprocessing of X-rays and if optical emis- 



sion originates in the disk itself due to its internal heat- 
ing. However, these two scenarios predict different lags 
between the optical and X-ray emission of IPs. In the for- 
mer case the optical emission should lag the X-rays, in the 
latter it should lead the X-rays. This issue can be studied 
using simultaneous optical and X-ray observations of the 
source. 

With a sample of magnetized white dwarfs we demon- 
strated that this behavior may be universal for white 
dwarfs with truncated accretion disks, like intermediate 
polars, with the break frequencies depending on the mass 
accretion rate. This behavior is very similar to that of ac- 
creting neutron stars with truncated accretion disks (e.g. 
accretion powered X-ray pulsars). 

We suggest that the values of the break frequencies in 
power spectra of accreting white dwarfs can be used to 
make estimates of the inner radii of the truncated accre- 
tion disks. In turn, these can be used to make estimates 
of the white dwarf magnetic fields. 

We proposed a few tests of the outlined paradigm. In 
particular, for observations of intermediate polars at dif- 
ferent mass accretion rates our interpretation predicts that 
the break frequency of their power spectra should increase 
with an increasing mass accretion rate in the binary sys- 
tem. 
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